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INTRODUCTION 


A  fundamental  requirement  for  the  calculation  of  transport 
and  diffusion  of  smoke  and  contaminants  in  the  I^wer  atmosphere 
is  the  specification  of  the  wind  field.  The  direct  approach  to 
computing  the  surface  winds  by  numerically  solving  the  equations 
of  motion  of  the  lower  atmosphere  is  impractical  due  to  both  com¬ 
putational  restrictions  and  lack  of  knowledge  of  many  of  the 
basic  processes  involved.  Most  currently  used  methods  for  com¬ 
puting  the  surface  wind  field  rely  upon  observations;  the  results 
produced  by  interpolating  between  observed  points.  The  horizon¬ 
tal  resolution  of  these  fields  is  limited  by  the  density  of 
observations . 

For  tactical  field  operations  there  is  a  need  for  a  model 
which  is  able  to  calculate  high  resolution  surface  wind  fields 
over  arbitrary  terrain  without  the  need  for  observations.  The 
model  should  account  for  the  specific  influences  of  day  of  the 
year,  time  of  day,  atmospheric  transmissivity,  surface 
temperature,  cloud  cover,  atmospheric  flow  and  surface  roughness. 
As  well,  it  should  be  computationally  efficient  so  that  it  can  be 
used  in  near  real  time. 

This  report  describes  the  development  of  a  surface  wind 
model  and  associated  terrain  processing  methods  designed  to  meet 
the  above  criteria.  The  method  is  an  extension  and  improvement 
of  a  method  originated  by  Ryan  (1974).  The  basic  principle  of 
Ryan's  work  is  to  empirically  incorporate  the  commonly  observed 
influences  that  terrain  has  on  surface  winds.  Second  order 
interaction  between  the  influences  is  ignored.  This  empirical 
approach  while  greatly  facilitating  the  computational  problem 
offers  realistic  winds  for  transport  and  diffusion  models.  In 
several  case  studies,  the  method  successfully  calculated  the  time 
and  space  variations  of  winds  in  complex  terrain. 

The  work  proceeded  in  two  primary  tasks.  Ryan  had  relied 
on  personal  measurement  of  terrain  for  use  in  his  method.  The 
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first  task  supplanted  this  approach  and  involved  developing  the 
methodology  necessary  to  produce  an  optimized  computation  grid 
from  high  resolution  uniformly  gridded  terrain  height  data.  The 
development  of  the  computation  grid  takes  into  account  the 
terrain  influences  on  the  wind  field  so  that  a  wind  field  of  any 
desired  resolution  can  be  obtained  for  all  terrain  types. 

The  second  task  involved  the  development  of  the  surface 
wind  model  which  incorporates  as  many  as  possible  of  influences 
terrain  has  on  the  local  wind.  The  model  calculates  wind  vectors 
on  the  optimized  (and,  normally,  nonuniform)  computation  grid 
produced  by  the  methodology  of  the  first  task.  The  resulting 
wind  field  can  then  be  interpolated  to  a  uniform  grid  of 
arbitrary  resolution  as  specified  by  the  user. 

The  details  of  the  algorithms  and  computer  programs  for¬ 
mulated  for  the  development  of  these  programs  is  discussed  in  the 
following  sections. 


SECTION  1 

TERRAIN  ANALYSIS  AND  GRID  GENERATION  (TAGG) 

The  objective  of  the  terrain  analysis  and  grid  generation 
computer  programs  is  to  provide  an  optimum  grid  and  associated 
data  for  the  computation  of  surface  wind  fields  for  a  given 
geographic  region.  The  system  analyzes  terrain  height  data  to 
select  appropriate  grid  points  and  then  computes  the  values  of 
the  terrain  wind  parameters  (TWPs)  for  these  points.  The  TWPs 
are  certain  data  items  -  slope,  aspect,  valley  direction,  etc. 
-required  for  calculating  local  influences  on  the  surface  wind. 
The  output  of  this  system,  the  optimum  {and,  most  likely, 
nonuniform)  grid  and  associated  TWPs,  forms  a  data  base  for  all 
surface  wind  computations  for  the  region  under  consideration. 

1.1  GENERAL  STRUCTURE 

The  modular  design  of  the  TAGG  system  is  shown  in  Figure  1 
The  components  are:  (1)  selection  of  the  analysis  region,  (2) 
Tchebyshev  polynomial  smoothing  of  the  terrain  height  field,  (3) 
calculation  of  the  terrain  wind  parameters  and  (4)  determination 
of  the  optimum  computation  grid. 

1.2  SELECTION  OF  THE  ANALYSIS  REGION 

This  initial  part  of  the  procedure  requires  user 
decision-making.  The  user  divides  the  large-scale  map  area  into 
overlapping  square  subregions  or  "terrain  windows."  Each  subre¬ 
gion,  'n,'  will  be  specified  by  a  length,  1.  An  overlap  between 
adjacent  windows  will  be  specified  by  the  parameter  al .  The 
overlap  is  required  for  consistency  in  the  smoothing  process 
which  follows.  These  parameters  are  constant  for  all  n=N  terrai 
windows.  The  UTM  coordinates  for  the  corners  of  the  large-scale 
map  will  also  be  given  by  the  user. 


Function 


Outpuc/Input 


Ogeraticm 


Subjective  Map 
Analysis 


T 


Number  of  windows:  (N) 

UTN  Corner  Coordinaeas:  (*c,yc) 
Window  Parameters:  (l,ai) 


User  selects  region  for 
analysis  and  size  of  sub¬ 
divisions  for  smoothing 


Q. 

Terri 


retrain  Tape 
Height  field  on 
uniform  fine  grid 


Output/ Input 


Output/ Input 


Wind  Error 
Tolerances 
‘min' ‘mas 


Data  Preprocessor 


Terrain  Heights  all 
arid  points.1  (hij>n 
for  Window,  n 


Tchebyshev  Smoother 
{hi3  !n*>(»t)n 


Expansion  Coefficients:  (at;n 
arid  Increment:  <4{,4it)n 
Fit  Efficiency:  fn 


TWP/Grid  Point 

Selector  Module 

Terrain  Wind  parameters:  (TWP)n| 
and  Grid  Points:  {*t.yi>n 


Selects  heights  in  each 
nth  Subdivision  (window) 


Smooths  height  field  in 
each  window  using 
Tchebyshev  polynomial 
series.  Tt(  ; ,  n) 


Computes  Terrain  Wind 
Parameters  and  selects 
optimum  nonuniform 
grid  points 


Output 


n>H 


n  *  n*l 


Figure  1.  General  Software  Structure  of  the  Terrain 
Gridding  System. 
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1.3  TCHEBYSHEV  POLYNOMIAL  SMOOTHING  OF  THE  TERRAIN 


The  fourth-order  two-dimensional  Tchebyshev  polynomial  was 
selected  to  fit  a  gridded  field  of  terrain  heights.  A  fourth 
order  polynomial  in  two  dimensions  may  have  up  to  three  critical 
values  (slope  =  0)  along  any  line  of  sight  through  the  repre¬ 
sented  region.  Thus  the  degrees  of  freedom  of  the  polynomial 
allows  for  two  peaks  and  a  valley  or  two  valleys  and  a  peak  to 
occur  in  the  smoothed  data  in  any  direction.  The  subjective 
choice  of  window  size  should  be  based  on  the  horizontal  dimen¬ 
sions  of  important  terrain  features.  No  more  than  two  meteorolo¬ 
gically  significant  hills  or  valleys  should  be  included  within 
the  dimension  of  a  single  window.  This  selection  enables  a  large 
number  of  terrain  heights  to  be  represented  by  a  single 
polynomial  fit.  In  addition,  the  fitting  removes  the  small  scale 
terrain  "bumps"  which  are  too  small  to  significantly  affect  the 
local  wind. 

Tchebyshev  polynomials  are  advantageous  for  at  least  two 
reasons.  First,  they  are  orthonormal,  thus  their  coefficients 
are  readily  obtained.  Secondly,  they  are  efficient  in  repre¬ 
senting  complex  terrain  because  they  have  infinite  order  con¬ 
vergence  for  non-periodic  data  fields  (Boyd,  1978).  The  latter 
property  is  the  chief  advantage  Tchebyshev  polynomials  have  over 
Fourier  series. 

The  smoothed  height  field  derived  by  fitting  fourth  order 
two-dimensional  Tchebyshev  polynomials  to  raw  terrain  height 
data  is 


h(x ,y) 


-  I  l 

m  n 


a 

m,n 


T 


m,  n 


where  the  Tchebyshev  polynomials 


T 

m,n 


T  (x) 
m 


Tn<y> 


(1) 


(2) 
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r 


are 


T0,0  *  1 
*1,0  *  x 
t2,0  *  2x2~1 
t3,0  “  4x3-3x 
t4,0  “  8x4“8x2+1 

*o’,i  =  y 
*1,1  =  xy 

*2,1  =  (2x2-l)y 
t3,1  =  (4x3-3x)y 

*0,2  =  2y2-i 
*i,2  ■  x(2y2-i) 

*2,2  -  (2x2-l)(2y2-l) 
T0, 3  =  <  4y3-3y) 

*1,3  =  x(4y3-3y) 

*0,4  -  8y4-8y2+l 


and  the  polynomial  coefficients  are  derived  from 


Jo  Jo  *u'y> 

°W,n  *  1-1  J-l  _ 


l  l  *, 

i=0  j=0 


m,n(xi,yj) 


where 

Z(x,y):  raw  terrain  height 

i:  index  in  x  direction 
j:  index  in  y  direction 

I,J:  total  number  of  data  points  in  x  and  y  directions 
respectively. 


After  the  smoothing  is  performed, the  residual  or  root  mean 
square  height  differences,  e,  between  the  raw  and  smoothed  height 
field  is  calculated.  The  residual  is  a  measure  of  the  fidelity 
with  which  the  polynomial  series  represents  the  original  height 
field.  An  efficiency  of  fit  for  each  window  is  obtained  from  the 
root  mean  square,  e,  and  the  standard  deviation  of  the  raw  height 
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field,  o 


as 


z' 


Qualitatively,  a  fit  efficiency  of  0.60  and  above  represents  an 
adequate  fit  while  an  efficiency  of  0.85  and  above  means  an 
excellent  fit. 

Subjective  map  analysis  is  recommended  so  that  the  window 
size  is  chosen  to  best  characterize  the  terrain  features.  The 
smoothing  which  occurs  as  a  consequence  of  the  fit  can  then  work 
to  the  user's  advantage  for  filtering  out  irrelevent  details  in 
the  terrain. 

Figures  2  and  3  compare  actual  and  fitted  height  fields  for 
two  idealized  terrain  surfaces.  Respectively,  they  are  a  cosine 
shaped  knoll  and  a  cosine  shaped  valley.  The  general  charac¬ 
teristics  of  each  surface  is  reasonably  conserved.  Note  that  in 
the  cosine  knoll  case,  which  is  surrounded  by  a  flat  region,  an 
upward  sloping  region  occurs  near  the  corners.  Figures  4  and  5 
compare  an  actual  terrain  height  field  (shown  as  the  deviation 
from  the  mean  elevation)  with  its  Tchebyshev  polynomial 
representation.  While  smoothing  is  evident,  the  general 
complexity  of  the  original  terrain  is  preserved  by  the  fitting. 
The  efficiency  of  fit,  fn,  for  this  case  is  0.822. 

In  practice  a  sizeable  number  of  individual  terrain  fits 
may  be  required  to  represent  an  area  of  interest.  A  mosaic  of 
overlapping  fit  windows  is  used  to  produce  a  relatively  con¬ 
tinuous,  smoothed  terrain  surface.  Centered  within  each  window 
is  a  region  of  computational  validity  where  the  wind  vectors  are 
calculated  called  a  cell  (see  figure  below).  The  window 
overlapping  prevents  sharp  discontinuities  that  might  otherwise 
occur  at  the  boundaries  between  adjacent  windows.  In  the  cosine 
knoll  case  for  example,  the  upward  sloping  regions  near  the 
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corners  would  not  be  used  in  wind  computations  since  they  would  not 
be  within  the  cell. 


Window  2 


I - 1 


Window  1 


1.4  TERRAIN  WIND  PARAMETER  (TWP)  GENERATION 

The  wind  field  at  each  grid  point  is  influenced  by  certain 
terrain  parameters  of  its  surroundings.  These  Terrain  Wind 
Parameters  (TWPs)  are  derived  from  the  Tchebyshev  polynomial 
representation  of  the  terrain  surface.  For  the  polynomial  repre¬ 
sentation  of  a  window,  TWPs  are  calculated  only  at  grid  points 
within  the  cell  boundaries  as  discussed  in  the  previous  section. 
The  TWPs,  since  they  describe  terrain  features,  need  to  be  calcu¬ 
lated  only  once  for  any  given  cell.  The  ability  to  calculate 
many  wind  fields  without  the  need  to  handle  the  basic  terrain 
height  data  each  time  greatly  reduces  computer  costs.  As  well, 
the  analytical  representation  of  the  terrain  surface  allows  for  a 
more  straightforward  determination  of  the  TWPs.  The  following 
paragraphs  discuss  the  definition  and  calculation  of  each  TWP. 

Local  Slope,  S. 

The  local  slope  S  at  a  grid  point  is  the  modulus  of  the 
local  terrain  gradient,  i.e.. 
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s  •  I  vh  i a  u  15  )2  +  C  *7  )V/2  (5) 

It  is  calculated  directly  from  the  polynomial  expression  for  the 
terrain  height  surface. 


Aspect  Angle,  a. 


The  aspect  angle  at  a  point  is  the  direction  towards  which 
the  tangent  plane  is  facing.  It  is  defined  such  that  a  slope 
(Plane)  facing  southward  has  an  aspect  angle  of  180°.  A  slope 
facing  northward:  0°.  The  aspect  angle  can  be  seen  to  be 
equivalent  to  the  compass  direction  opposite  to  the  direction  of 
the  terrain  gradient  vector,  7h.  Hence 


a 


ir/2  -  tan 


( 


3h/ty  , 
3h/3x  >  ' 


if  a  <  0  a  =  a  +  2* 


(6) 


Upwind  and  Downwind  Slope  to  Horizon,  Y. 

At  each  grid  point  (x,y)  with  smoothed  height  coordinate, 
h,  the  slope  to  horizon,  Y,  is  calculated  in  the  eight  primary 
compass  directions.  We  use 


Y  =  Uh/As)max 


(7) 


where  Ah  =  h  -h 
P 

and  As  =  [ ( x 

P 


X)2  +  (yp  -  y)2]^2 


(8) 

(9) 


The  slope  is  calculated  between  the  point  (x,y)  and  every  other 
grid  point  (Xp,  yp)  in  the  specified  direction,  the  maximum  value 
of  Ah/ As  along  a  specified  compass  direction  defines  the  horizon. 
The  eight  compass  directions  are  defined  to  be  along  lines  of 
constant  x,  constant  y,  and  x=±y.  The  coordinate  system  is 
defined  such  that  increasing  x  is  eastward  and  increasing  y  is 
northward.  Depending  on  the  direction  of  the  general  wind,  the 
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appropriate  upwind  and  downwind  values,  and  ,  are  selected 
from  the  eight  values  in  the  set  {Y>. 

Ridge  Parameters,  Dh,  Ds,  Hr. 

Ridge  parameters  are  calculated  for  those  points  that  lie 
in  the  vicinity  of  a  ridge.  More  explicitly,  the  parameters  Dh, 
Ds,  and  HR  are  calculated  for  those  points  that  lie  between  the 
zero  curvature  point  (concavity  change  from  valley  to  ridge)  and 
the  ridge  line  point.  The  ridge  line  point  in  a  specified  direc- 
"  tion  is  derived  from  the  gradient  of  the  smoothed  terrain  eleva¬ 
tion  in  that  direction.  The  gradient  is  a  cubic  polynomial  in 
x,y  coordinates  and  is  then  transformed  into  polar  r,  #,  coor¬ 
dinates  with  origin  at  the  grid  point  of  calculation.  Thus 


r*Vh  =  cos  <j>  ^  +  sin  $  ay  =  Ao  +  Alr  +  A2f2  +  A3r3  (10) 

where  r  =  i  cos  $  +  j  sin  $  :  unit  radial  vector 
<1>  :  clockwise  angle  from  North 
r  :  radial  coordinate 

The  real  roots  of  the  above  expression  are  evaluated  for  each  of 
eight  compass  directions. 

The  parameter  Dh  is  the  elevation  difference  between  the 
zero  curvature  point,  hz,  and  the  ridge  line  in  a  specified 
direction,  that  is: 


Dh  ■  |  hR  -  h2  |  Ul) 

where  hr  :  smoothed  value  of  elevation  at  the  ridge  line. 

The  parameter  Ds  is  the  straight  line  horizontal  distance 
from  the  calculation  point  to  the  ridge  line  in  a  specified 
direction,  i.e. 

s  -  [(x  -  xR)2  +  (y  -  yR)2]^  (12) 
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where  xR,yR  :  coordinates  of  the  point  on  ridge  line  in  the 
specified  direction. 

The  roots  {r^}  of  the  cubic  correspond  to  either  a  point  on 
the  valley  line  or  ridge  line  in  the  direction  $  from  the  grid 
point  of  calculation.  A  valley  line  point  is  discriminated  from 
a  ridge  line  point  by  the  sign  of  the  partial  derivative  of  the 
gradient. 


r  •  Vh) 


(Aj^  + 


2A2r 


3A3r2) 


r=r 


<  0,  Ridge 
>  0f  Valley 


(13) 


The  parameter  HR  is  a  measure  of  the  ridge  width  in  a  spec¬ 
ified  direction.  The  ridge  width  is  the  distance  in  the  hori¬ 
zontal  plane  between  the  ridge  line  point  and  the  zero  curvature 
point  (x  ,  y  ) .  The  zero  curvature  point  is  the  real  root  of  the 
equation, 

A1  +  2A2C  +  3A3r2  =0  (14) 

that  lies  between  the  calculation  point  and  the  valley  line 
point.  The  ridge  width  is  calculated  in  each  of  up  to  eight 
directions  as: 

Hr  =  t(xR  ’  V2  +  (yR  "  yz)2]^  <15) 

The  ridge  parameters  are  used  to  calculate  enhanced  wind 
speed  for  grid  points  lying  on  the  upper  half  of  a  hill.  No 
ridge  parameter  are  calculated  for  points  on  the  lower  half 
(i.e.,  below  the  zero  curvature  point  in  a  valley).  Ridge  points 
are  calculated  in  each  of  eight  compass  directions  when  possible. 
In  many  cases  a  ridge  point  will  not  be  found  in  certain  direc¬ 
tions  in  which  case  no  parameters  are  calculated. 

General  Valley  Direction,  0. 

The  general  valley  direction  0,  is  the  upvalley  direction 
of  the  "dominant"  valley  within  the  cell  (non-overlapping  part) 
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of  a  window.  The  choice  of  resolution  or  window  size  used  in 
fitting  the  Tche.^yshev  polynomials  to  the  original  grid  of 
terrain  heights  should  be  such  that  the  cell  size  is  on  the  order 
of  the  dimensions  of  a  valley.  At  each  grid  point  within  the 
cell  the  average  valley  direction  of  the  entire  cell  is  assigned. 
The  average  direction  is  the  vector  mean  of  the  individual  valley 
directions  at  each  grid  point.  The  calculation  of  an  individual 
valley  direction  at  a  point  begins  by  looking  in  the  four  direc¬ 
tions  (N,E,S,W)  for  the  closest  valley  point.  The  valley  points 
are  determined  from  the  roots  of  the  polynomial  in  a  similar 
fashion  to  finding  a  ridge  point.  At  the  closest  valley  point  a 
search  is  made  in  10°  increments  for  the  direction  of  greatest 
elevation  decrease.  The  individual  valley  direction  at  the  grid 
point  is  then  defined  to  be  180°  from  the  direction  of  greatest 
elevation  decrease. 

Valley  Concavity,  Y. 

The  valley  concavity,  Y  is  a  measure  of  the  "valleyness"  of 
the  valley  at  an  arbitrary  grid  point.  The  equation  for  valley 
concavity  is 

Y  =  <Yi  +  Yr)/2 

where  Y^  and  Y^  are  the  slopes  to  the  horizon  to  the  left  and 
right  of  the  general  valley  direction.  The  slopes  to  the  horizon 
are  derived  the  same  way  as  the  general  slope  to  the  horizon 
( Equation  7) . 

1.5  DETERMINATION  OF  THE  OPTIMUM  WIND  COMPUTATION  GRID 

The  construction  of  an  optimum  computation  grid  is  a  two 
step  process.  The  first  step  is  the  construction  of  a  uniform 
grid  of  sufficient  resolution.  The  second  step  is  to  remove  all 
unnecessary  grid  points  from  the  uniform  grid.  The  resulting 
grid  is  nonuniform,  due  to  the  removal  of  some  of  the  uniform 
grid  points,  and  optimal  since  no  more  points  can  be  removed. 

The  overall  process  is  driven  by  three  user  specified  input 
values  emax'  emin'  an(3  Ne* 
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In  the  first  step  of  constructing  a  uniform  grid  of  suf¬ 
ficiently  fine  resolution,  the  two  values  emax  and  Ne  are  used. 
emax  is  nominal  maximum  permissible  value  of  wind  speed  dif¬ 

ference,  AV,  between  any  two  adjacent  grid  points.  However,  this 
restriction  is  relaxed  to  the  extent  that  up  to  Nc  values  of  AV 
may  exceed  enax. 

The  uniform  grid  is  constructed  in  the  following  manner. 

The  program  begins  by  constructing  a  four  by  four  grid  of  test 
points  to  represent  a  terrain  cell  (  a  cell  is  the  nonoverlapping 
portion  of  the  Tchebyshev  window).  The  program  then  computes  the 
terrain  wind  parameters  at  each  of  the  sixteen  grid  points.  For 
each  grid  point  the  program  then  computes  the  value  of  AV  (see 
Section  2)  between  it  and  its  right  and  top  neighbor  (except  for 
the  right  column  and  top  row  of  points).  If  a  value  of  AV 
exceeds  emax  a  counter  is  incremented  and  a  check  is  made  to  see 
if  the  counter  values  exceeds  Ne.  If  this  happens,  then  the 
program  "knows"  the  grid  resolution  is  not  fine  enough  (i.e.,  the 
number  of  AV,  exceeding  emax  >  N£).  In  such  a  case  the  program 
doubles  the  grid  resolution  (four  by  four  becomes  an  eight  by 
eight)  and  starts  the  process  again. 

If  after  computing  AVs  for  all  the  grid  points,  the  number 
of  "too  big"  AVs  does  not  exceed  Ne  then  the  program  accepts  the 
current  grid  resolution  as  being  sufficiently  fine.  It  should  be 
noted  that  the  program  will  continue  doubling  the  grid 
resolution,  if  necessary,  up  to  a  maximum  of  64  by  64  grid 
points.  If  such  a  resolution  is  still  not  fine  enough,  the 
program  will  print  an  error  message  proclaiming  this  condition 
and  continue  on  with  the  grid  optimization  procedure.  The 
program  will  utilize  the  64  by  64  grid  even  though  this  grid  does 
not  satisfy  the  emax  and  Ne  criteria. 

Once  a  fine  enough  uniform  grid  has  been  constructed,  the 
program  proceeds  to  remove  all  unnecessary  grid  points.  This 
process  is  controlled  by  the  user  input  value  Em^n.  A  grid  point 
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is  deleted  from  the  grid  if  the  values  of  AV  between  the  point 
and  its  upper#  right,  lower,  and  left  neighbors  are  less  than 
emin*  Each  point  in  the  grid  is  tested  for  this  condition  and 
deleted  if  possible.  After  the  first  deletion  pass  through  the 
grid,  some  of  the  remaining  grid  points  have  new  neighbors  due  to 

the  deletion  of  their  original  neighbors.  A  second  deletion  pass 

% 

is  made  through  the  now  sparser  grid  to  delete  any  of  the 
remaining  grid  points  that  also  satisfy  the  deletion  criteria 
with  their  new  neighbors.  Deletion  passes  through  the  grid  are 
repeated  until  a  deletion  pass  is  made  in  which  no  grid  points 
are  deleted.  This  means  that  no  grid  points  get  new  neighbors, 
thus  no  new  aVs  will  be  computed  which  could  be  less  than 
emin  and  ths  deletion  process  is  completed.  The  grid  points  that 
remain  comprise  the  optimal  nonuniform  grid. 

It  should  be  noted  that  the  deletion  process  can  be  "turned 
off"  by  the  user  by  entering  an  input  value  of  0  for  emin*  Since 
aV  is  always  greater  than  or  equal  to  0,  no  points  will  be 
deleted.  Likewise  the  user  can  "force"  the  program  to  accept  a 
four  by  four  grid  every  time  by  entering  a  value  Ne  greater  than 
32.  Since  the  program  only  checks  the  values  of  AV  to  the  top 
and  right  neighbors  of  a  grid  point,  a  maximum  of  32  "too  big" 

AVs  can  occur. 

To  test  and  evaluate  the  grid  optimization  procedure,  a 
synthetic  terrain  data  base  was  generated  that  represents  a  three 
dimensional  cosine  curve.  To  test  the  optimization  procedure  the 
wind  difference  criteria,  AV,  was  replaced  by  the  easily  calcu¬ 
lable  height  difference  between  the  points.  Because  this  was  a 
test  of  the  optimization  process  and  not  of  the  grid  construction 
process,  a  grid  of  130  by  130  points  was  constructed,  which 
correspond  to  an  emax  of  .2  altitude  units.  In  this  test  case, 
the  program  increased  the  qrid  resolution  by  ten  times  at  each 
step  rather  than  doubling  the  resolution.  Figures  6,  7,  and  8 
show  this  grid  after  optimization  corresponding  to  cm^n  values 
of  .05,  .1,  and  .2  respectively.  It  can  be  observed  from  these 
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figures  that  an  increasing  value  of  em^n  will  cause  the  resultant 
optimal  grid  to  thin  out.  Note  that  the  areas  with  no  points 
correspond  to  the  peaks  and  troughs  of  the  cosine  curve  where  the 
altitude  difference  between  points  is  small. 
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SECTION  2 


THE  SURFACE  WIND  MODEL 

The  Surface  Wind  Model  (SWM)  calculates  the  surface  winds 
over  a  specified  region  on  the  optimal  computation  grid.  In 
order  to  make  wind  calculations  at  each  grid  point  it  is 
necessary  to  know  the  terrain  wind  parameters  (TWPs)  and  surface 
roughness  at  each  grid  point,  the  latitude,  longitude,  cloud 
cover  and  transmissivity  for  the  specified  region,  and  the  time 
and  date.  The  procedure  for  using  these  parameters  to  derive 
surface  winds  is  described  below.  The  resulting  wind  field  on 
the  nonuniform  grid  is  interpolated  to  a  uniform  grid  for  final 
presentation. 

The  SWM  consists  of  three  main  components  (Figure  9).  The 
first  component  involves  the  calculation  of  the  atmospheric  sta¬ 
bility  and  general  surface  wind  from  the  free  atmospheric  wind 
(at  about  the  1500  m  or  the  850  mb  level).  The  general  surface  wind 
is  obtained  for  each  square  kilometer  of  the  analysis  region. 

The  resulting  wind  vectors  are  passed  to  component  2  which  calcu¬ 
lates  the  six  terrain  perturbations  to  the  general  surface  wind. 

The  perturbations  are  calculated  at  each  grid  point  using  the 
TWPs  associated  with  that  point.  The  six  perturbations  are  vec- 
torially  added  to  the  general  surface  wind  for  the  appropriate 
square  kilometer  to  give  a  resultant  local  surface  wind  at  each 
grid  point.  The  third  component  of  the  surface  wind  model,  the 
interpolator  converts  the  wind  field  from  a  nonuniform  grid  to  a 
use  specified  uniform  grid  for  output  and  display  purposes. 

2.1  INPUT 

Input  to  the  surface  wind  model,  as  shown  in  Table  1,  comes 
from  three  sources:  the  terrain  and  TWP  output  from  the  terrain 
orid  analysis,  aerodynamic  surface  roughness  values  for  each 
square  Kilometer  region,  and  input  describing  the  meteorological 
condition  for  the  time  of  wind  calculations.  The  optimum  non- 
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Figure  9.  Surface  Wind  Model. 
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uniform  grid,  the  raw  terrain  height,  and  the  TWp's  for  each  grid 
point  are  direct  outputs  from  the  terrain  analysis  program.  The 
terrain  surface  roughness  value  for  each  square  kilometer  is 
obtained  independently  from  an  analysis  of  the  type  vegetation 
and  terrain  surface  for  the  region.  MacArthur  and  Haines  (1981) 
discuss  techniques  for  calculating  surface  roughness.  Other 
meteorological  input  necessary  to  run  the  surface  wind  model  is 
j  '  also  shown  in  Table  1.  This  input  includes  the  free  atmospheric 

I  wind,  temperatures,  heights  and  surface  temperatures  on  a  four 

j  point  grid  bounding  the  region  of  interest.  In  addition,  lati¬ 

tude,  time  of  day,  day  of  year,  percent  cloud  cover,  sea  breeze 
data  (if  appropriate),  atmosphere  transmissivity  and  program 
control  variables  are  also  required. 


The  following  sections  describe  the  three  main  components 
of  the  surface  wind  model:  the  calculation  of  the  general  sur¬ 
face  wind,  the  local  terrain  influences  on  the  wind,  and  the 
interpolation  of  the  wind  field  to  a  uniform  output  grid. 


TABLE  1 

INPUT  FOR  SURFACE  WIND  MODEL 


Input  Supplied  From  Terrain/Wind  Data  Base 


•  Wind  Computation  Grid,  (xi'yi'zi}  i  =  1,  N.  This 

is  the  nonuniform  grid  computed  by  the  Terrain  Gridding 
Package 

•  Terrain  Wind  Parameters  (TWP's)  at  each  point  i  on  the 
wind  computation  grid 

•  Surface  Roughness  Field.  Aerodynamic  surface  roughness 
length  for  each  square  km  of  the  analysis  region 

,  P  P  P 

lxi »  Yi>  z0i}'  i  =  1»  •••/  M  <  2500 


Input  Supplied  By  The  User 

•  Switch  to  ignore  initialization  and  allow  direct  input 
of  general  surface  wind  components  Ub>  v^. 

•  UTM  coordinates  of  the  four  corners  of  the  region  of 

interest  {x^,  yjJ,  i  *  1,  4.  Alternatively,  the 

SW  corner  point  (xc,  yc)  and  the  length  or  width  of  the 
region. 

•  Free  Atmosphere  (850  mb)  Wind,  Temperature,  and  Heights, 
and  surface  temperatures  on  a  four  point  grid  bounding 
the  region  of  interest 

a  a  a  a  a  a  a 
{ x i »  T^,  ZH^,  Tgfc^,  }  i  =»  1 ,  •«.,  4 

•  Latitude,  $r,  Time  of  day,  t<j ;  Day  of  Year,  Dy,  for 
analysis  region 

•  Percent  Cloud  Cover,  Cc,  and  Transmissivity,  PQ.  Used 
to  determine  stability. 

•  Sea  Breeze  Data.  Distance  from  analysis  region  (grid) 
center  to  shore  of  ocean  (lake,  etc.),  d^,  and  direction 
to  shore,  b^.  Direction  is  given  in  degrees  counter¬ 
clockwise  from  East. 

•  Program  Control  Variables.  Switches  to  activate  the 
various  modules  containing  the  correction  effects 
(sheltering,  diverting,  etc.) 
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Output  Grid.  Uniform  grid  for  presenting  the  output 
wind  field.  Given  as  the  SW  corner  coordinates  (xQ,  yQ), 
the  mesh  length  and  width  dimensions  (ax0,  Ay Q) ,  and  the 
overall  grid  dimensions  (Lq,  Wq).  Also  number  of  inter¬ 
polation  scans/  initial  radius  of  influence,  percent 
reduction  in  influence  radius  for  each  succeeding  scan, 
and  the  convergence  criterion,  which  when  satisfied 
terminates  the  interpolator. 


2.2  THE  GENERAL  WIND  AT  THE  SURFACE 


The  free  atmosphere  wind,  at  an  altitude  significantly 
above  the  influence  of  the  terrain  serves  as  the  basic  wind  input 
parameter  to  the  model.  From  the  free  atmosphere  wind  a  surface 
wind  field  at  6  meters  altitude  is  generated  for  each  square 
kilometer.  This  general  wind  at  the  surface  takes  into  account 

% 

influences  of  the  ground  in  flat  terrain.  Nonuniform  terrain 
influences  are  accounted  for  in  the  next  section. 

Variation  of  wind  speed  and  direction  between  ground  and 
the  free  atmosphere  is  dependent  upon  surface  characteristics 
(such  as  heating  or  cooling  and  roughness),  the  coriolis  force, 
and  the  temperature  structure  of  the  intervening  layer.  The 
following  scheme  combines  two  methods  to  derive  the  6  meter 
winds.  The  method  of  Gerrity,  1967,  is  used  to  deduce  the  50 
meter  wind  from  the  free  atmospheric  wind,  and  the  Goodin  and 
McRae  (1980)  method  is  used  for  the  wind  speed  profile  between 
the  ground  and  50  meters.  The  scheme  is  versatile  since  it 
incorporates  both  baroclinic  shear  as  well  as  the  effect  of  sur¬ 
face  roughness. 

Free  atmosphere  wind,  temperature,  height  and  surface 
temperature  are  specified  on  a  four  point  grid  surrounding  the 
analysis  region.  These  free  atmosphere  parameters  are  opera¬ 
tionally  available  from  National  Weather  Service  or  Air  Force 
Global  Weather  Center  weather  charts.  The  following  additional 
variables  are  needed  to  calculate  the  general  wind  at  the 
surface : 

•  cloud  cover  and  transmissivity  for  the  region  specified 
by  the  four  points 

•  surface  roughness  for  each  sq.  km  in  the  analysis 
region. 

•  latitude,  longitude,  time  of  day,  and  day  of  the  year. 

First,  the  geostrophic  wind  components  at  the  top  of  the 

constant  stress  layer  and  the  Ekman  turning  of  the  wind  between 
the  free  atmosphere  and  constant  stress  layer  is  obtained. 
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Gerrity's  model  divided  the  region  between  the  earth's  surface 
and  the  free  atmosphere  into  two  regions  as  shown  in  Figure  10. 
One  is  the  constant  stress  layer  near  the  earth's  surface  up  to 
50m  where  viscous  forces  dominate  and  the  other  is  a  transitional 
layer  between  the  constant  stress  regime  and  the  free  atmosphere. 
In  the  free  atmosphere  pressure  and  rotational  forces  usually 
balance.  Gerrity  proposed  that  the  geostrophic  wind  in'  the  tran¬ 
sitional  layer  could  be  modeled  as 


0  =  U 

g  9 

V  =  v  1 

g  g 


H 


+  B(H-z) 
+  C(H-z) 


(16) 

(17) 


where  U 


Vg  are  the  geostrophic  components  in  the  transition 


H  H 

layer  and  are  the  geostrophic  components  at  the  top  of 

the  transitional  layer  whose  height  is  H.  B  and  C  are  coef¬ 
ficients  evaluated  from  the  temperature  field  in  the  transitional 


layer  and  z  is  the  altitude, 
be  written  as 


More  explicitly  (16)  and  (17)  can 


„  .  u  “ .  u  -  (VM.  p.  f 

*  9  9  \  TH  /  *7  cJ 


3  1  a, 

3y  T  d 


(18) 


V  =  V 


+  V. 


9T, 


li 

f 


3  1  h  ■» 

3x  T  d 


(19) 


where  T  and  T„  are  the  temperatures  at  the  respective  heights, 

C  H 

of  50  meters,  and  H,  g  is  gravity,  and  f  is  the  coriolis  force 


for  the  latitude  of  interest. 


Th  is  simply  the  free  atmosphere 


temperature  while  Tc  is  50  meter  temperature  derived  from  the 
surface  temperature  and  an  appropriate  lapse  rate. 

A  linear  variation  of  temperature  in  the  transitional 
layer  is  assumed.  Thus,  the  vertical  temperature  profile  is. 


T ( z )  =  (T  -  T  ) ( z  -  c)  /  (H  -  c)  +  T 

rl  C  C 


(20) 
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Free  Atmosphere 


Figure  10.  Two  Part  Division  of  the  Planetary  Boundary  Layer. 


Substitution  of  (20)  into  the  integrals  in  (18)  and  (19)  yields 


<c  =  aJlnlVT=l  /  <th  -  V 


3T  3T 

+  lH  - cl  (T°iir  -  Tht*£)  >  Vo  (Th  -  Tc> 


3T  3T  3TU  3  T  _ 
-  <H  -  C)ln[TH/Tc]  /  l-af  *  (Jy-  -  Jy~> 


(21) 
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3T 
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3T 
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(22) 


Resubstituting  (21)  and  (22)  into  (18)  and  (19)  and  eva¬ 
luating  at  z  =  c  gives  the  geostrophic  wind  components  at  the  top 
of  the  constant  stress  layer. 

The  geostrophic  wind  at  the  top  of  the  constant  stress 
layer  can  be  corrected  for  the  sea  breeze  influence.  The  magni¬ 
tude  of  this  effect  depends  on  distance,  d.  ,  and  direction  b.  to 

b  b 

the  ocean.  Other  details  concerning  the  sea  breeze  remain  as  in 
the  original  Ryan  model  (1974). 

The  angular  deviation,  of  the  actual  wind  from  the 
geostrophic,  as  the  top  of  the  constant  stress  layer  is  based  on 
an  empirical  correlation  with  the  surface  Rossby  number,  R  , 
Gerrity  (1967) , 


\|>  3  a  (log  R  )  +  b  (log  R  )  +  c 

o  o 


(23) 


where 
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a  *  0.625, 

b  »  -12.75, 

c  *  80.625, 

R  =  G/f,z  ,  and 
o  o  o 

G  a  magnitude  of  geostrophic  wind  at  50m 

f  *  coriolis  force 
o 

zq  =  surface  roughness  height 


(U  2  +  V  2)J^. 

g  g 


Table  2  below  shows  the  variation  of  y  with  G  and  zQ  for 
f  =  10  /sec.  Equation  (23)  is  used  to  determine  the  wind 
direction  at  the  surface  from  the  direction  of  the  free 
atmosphere  wind. 


Equation  (23)  is  used  to  correct  the  wind  direction  at  the 
surface  (since  no  additional  turning  is  considered  to  occur  in 
the  constant  stress  layer)  for  the  influence  of  friction.  It  is 
currently  formulated  only  for  a  neutral  boundary  layer.  Angular 
deviation  would  be  larger  or  smaller  for  respectively  convective 
or  stable  boundary  layers. 


TABLE  2 


Angular 

Deviation 

of  the  Constant 

Stress 

Layer  Wind  from  Geostroph 

*# 

Ro 

G(m/sec) 

z0(m) 

Description 

39.6° 

104 

1 

1 

rough,  light  winds 

34.5° 

5 

x  104 

5 

1 

rough,  moderate  winds 

32.5° 

105 

10 

1 

rough,  windy 

26.6° 

106 

1 

10-2 

smooth,  light  winds 

23.3° 

5 

x  106 

5 

10"2 

smooth,  moderate  winds 

22° 

107 

10 

10"2 

smooth,  windy 

The  wind  profile  magnitude  within  the  constant  stress  layer 
is  based  on  stability  and  roughness.  A  quantitative  measure  of 
the  stability  influence  (Liu  et  al.,  1976)  is  based  on  the  stress 
layer  wind  and  insolation.  Briefly  a  is  calculated  as 

o  *  -1/2(3  -  Cw  +  |  Ce  |  )  Sign  (Ce)  (24) 


32 


where 


Cw  * 


Sign  (Ce)  = 


Ce  = 


G/2, 

4, 


0  £  G  £  8m/sec 
G  >  8m/sec 


(25) 


1,  Ce  >  0 

0,  Ce  =  0  (26) 

-1,  Ce  <  0 

3  strong 

2  moderate  day  time  insolation 

1  slight 

(27) 

0  heavy  overcast  }  day  or  night 
-1  4/8  cloud  cover 

night  time 

-2  <  3/8  cloud  cover 


The  algorithm  for  daytime  Ce  values  is  based  on  solar  angle 
and  cloud  cover.  The  representation  for  o  is  related  to  the 
Monin-Obuhkov  length,  L,  and  surface  roughness  zq,  by  Liu  et.  al . , 
as 


where 


L 


b2  M  ♦  b3”2) 
al°  +  a20^ 


(28) 


ax  =  0.004349 
a2  =  0.004724 

If  o  *  0r  then  L  *  106m. 


bx  -  0.5034 
b2  =  0.2310 
b3  =  0.0325 
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Following  Goodin  and  McRae  (1980),  the  surface  friction 
velocity,  u* ,  is  derived  from  the  geostrophic  wind,  G,  at  the 
reference  height  zr  (50  meters)  as 


* 

u 


(29) 


where  k  is  Von  Kantian' s  constant  (0.35)  and 


where  <)>(z/L)  a  function  related  to  wind  shear  is 
„  ,  -1/4 

♦  !  *  u  -  is(  1)) 

The  integral  form  of  the  velocity  gradient,  Goodin  and  McRae 
(1980)  is 


u(z)  *  u(zr)  + 


u*  T 

K  J 


rz  \  az 

'“L J  Z 


which  integrated  is 
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Using  z  =  6  meters,  z^  *  50  meters,  we  derive  our  general 
surface  wind  speed  from  equation  31.  Derivation  of  the  general 
surface  wind  is  summarized  in  Figure  11. 


2.3  LOCAL  INFLUENCES  ON  THE  GENERAL  WIND 

Predicting  the  surface  wind  speed  and  direction  in  complex 
terrain  is  a  particularly  difficult  problem.  The  resultant  wind 
at  a  given  location  arises  from  the  interaction  of  different 
scales  of  motion-synoptic,  mesoscale,  and  local.  The  problem  can 
become  quickly  intractable  if  the  interactions  between  various 
scales  of  motion  are  considered.  To  render  the  problem 
tractable,  the  assumption  is  made  that  second  order  interaction 
between  various  scales  of  motion  can  be  neglected. 

Diverse  mechanisms  are  at  work  influencing  the  wind  com¬ 
ponents  in  complex  terrain.  Synoptic  scale  influences  which  have 
already  been  considered  arise  from  the  large  scale  pressure 
gradient,  the  coriolis  force,  and  surface  friction.  A  mesoscale 
influence,  the  sea  breeze,  arises  from  land  and  sea  temperature 
differences.  The  sea  breeze  is  influenced  by  coriolis  force  and 
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Figure  11.  Steps  in  Calculating  General  Surface  Wind  from  Free 
Atmosphere  Kind. 
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friction  as  well.  Nonuniform  terrain  features  effect  both  synop¬ 
tic  and  mesoscale  flows  by  sheltering  and  diverting.  The  mean 
flow  sheltering  and  diverting  consist  respectively  of  a  decrease 
in  wind  speed  in  the  lee  of  obstacles  and  a  turning  of  the  wind 
when  it  is  obstructed.  Synoptic  scale  and  mesoscale  flows  are 
also  affected  in  other  ways  by  terrain.  The  flow  over  a  ridge  or 
hill  speeds  up  according  to  the  obstacle  shape  and  slope.  Flow 
separation  over  both  ridges  and  hills  is  dependent  on  the  terrain 
slope  and  atmospheric  stability,  but  occurs  at  different  values 
for  hills  (three-dimensional)  and  ridges  (two-dimensional).  The 
extent  of  separation  beyond  the  obstacle  is  also  different  for 
hills  and  ridges.  Nonuniform  terrain  also  exerts  local  influen¬ 
ces  on  the  wind  field.  Slope  and  valley  winds  arise  from  the 
temperature  differences  due  to  uneven  heating  or  cooling  between 
the  ground  and  air  away  from  valley  walls  or  mountain  slopes. 

Wind  component  variations  over  nonuniform  terrain  also 
occur  in  time.  Diurnal  and  seasonal  variations  must  be 
considered.  Thus,  the  resultant  wind  in  complex  terrain  consists 
of  many  terrain  induced  wind  perturbations  whose  amplitudes  vary 
in  both  time  and  space.  A  model  for  the  local  wind  must  incor¬ 
porate  both  the  space  and  time  variations.  The  following  sec¬ 
tions  describe  each  terrain  induced  perturbation  and  its 
influence  at  the  6  meter  level. 

2.3.1  Sheltering  and  Diverting 

The  general  wind  flow  at  a  point  in  complex  terrain 
is  modified  by  sheltering  and  diverting.  The  decrease  of  wind 
speed  on  the  lee  side  of  a  ridge  at  the  6  meter  level  due  to 
sheltering  is  computed  by  a  sheltering  factor,  Fu,  (Ryan,  1974): 

Fy  »  tan-1( .17YU)/100  .  (32) 

where  Yy  is  the  upwind  slope  to  the  horizon.  This  factor  is 
based  on  Van  Eimern's  (1955)  study  of  wind  reduction  by  the 
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Harburg  and  Geest  hills  near  Hamburg  and  Kaiser's  (19S9)  study  of 
wind  speed  reduction  by  shelterbelts. 


The  accompanying  figure  (12)  lends  additional  sup¬ 
port  to  Equation  (32).  Shown  are  the  Van  Eimern  results  for 
reduction  of  the  wind  speed  in  the  Elbe  Valley  near  Hamburg  and 
the  computed  factor  Fu  from  Equation  (32).  A  curve  was  fit  to 
this  data  that  nearly  parallels  the  Fu  curve.  Also  shown  are 
Counihan's  (1969)  results  for  the  turbulent  wake  behind  a  rect¬ 
angular  block,  which  also  favorably  compares  to  Fa. 


by 


The  surface  wind  speed  is  modified  for  sheltering 


W 


sd 


*  wb  -  lwb 


V 


(  33) 


where  is  the  sheltered  wind  speed  and  is  the  general  sur¬ 

face  wind  speed  from  equation  (31). 

The  wind  direction  is  modified-  by  a  diverting 

angle,  F^,  which  depends  on  the  surface  wind  direction  at  6 

meters,  0.  ,  aspect  a,  and  slope  Y.  of  the  terrain  downwind, 
o  a 

From  Ryan  (1974),  the  equation  is: 


Fd  =  -22.5Ydsin[2(a  -  o  )].  (34) 

Diversion  is  greatest,  22.5°,  for  a  slope  of  100%  and 
an  angle  between  aspect  and  wind  direction  of  45°.  The  diversion 
angle  goes  to  zero  as  the  wind  either  blows  along  or  normal  to 
the  slope.  The  diverting  angle  corrects  the  general  wind  direc¬ 
tion,  eb,  by 


sd 


-  eb  +  Fd 


(35) 


Following  Ryan,  the  east-west,  U  and  north-south,  V  , ,  com- 

sd  sd 

ponents  of  the  surface  wind  modified  for  sheltering  are 


sin  0  . 
sd 

COS  0  , 
sd 


(36) 

(37) 
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Van  Eimern  measurements 


SLOPE  TO  HORIZON  UPWIND  IN  PERCENT 


We  have  associated  a  relative  confidence  factor 
with  each  local  correction  that  is  applied  to  the  general  surface 
wind.  This  confidence  factor  indicates  the  degree  to  which  we 
believe  the  semi-empirical  equations  explain  the  actual 
phenomena.  The  relative  confidence  in  the  sheltering  and 
diverting  components  of  the  wind  is  quite  high,  thus,  we  have 
weighted  Ugd  and  Vgd  by  1.0. 

2.3.2  Slope  Wind 

In  mountainous  terrain  either  upslope  or  downslope 
winds  are  frequently  observed.  These  winds  arise  from  unequal 
heating  or  cooling  of  slope  surfaces  in  comparison  to  air  away 
from  the  slope  at  an  equal  altitude.  The  heating  promotes  posi¬ 
tive  buoyancy  by  day  and  negative  buoyancy  by  night.  Other 
influences  such  as  adiabatic  changes,  slope  of  terrain,  and 
nature  of  the  slope  surface  affect  the  strength  of  this  flow. 

According  to  Petkovsek  and  Hocevar  (1971),  the  night¬ 
time  slope  wind  strength  parallel  to  a  slope  is 


W 


ps 


( Ya  -  Y'J  (sin(s)) 


U  ~  exPifr'(Y< 


-  y')sin  SAt]  } 


where 


=  slope  wind  speed  parallel  to  slope, 

=  diabatic  cooling, 

=  adiabatic  lapse  rate, 

=  atmospheric  lapse  rate  outside  cooled  layer, 
=  slope  to  horizon 
=  acceleration  of  gravity, 

=  surface  friction  coefficient, 

=  temperature  outside  cooled  layer, 

*  time  since  sunset  or  sunrise. 


(38) 
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Petkovsek  and  Hocevar  claim  that  winds  calculated  from  (38)  match 
data  from  observations  well.  Petkovsek  and  Hocevar1 s  term 
"(y  -  y')“  better  models  actual  atmospheric  lapse  rate  con- 

a 

ditions  than  that  of  an  earlier  version  (Fleagle,  1950)  which 
considered  the  atmosphere  outside  the  slope  wind  layer  to  be 
isothermal.  Such  a  condition  restricts  the  maximum  values  pre¬ 
dicted  by  (38). 

Ryan  extended  (38)  to  the  daytime  upslope  condition 
by  making  the  term  apply  to  surface  heating  or  cooling  for 
respectively  day  and  night.  Ryan  considered  only  the  component 
of  the  slope  wind  in  the  horizontal  plane.  His  equation  is: 

Ws  =  2(1  +  Cn)  (1  -  Cp)  /  tanS,  (39) 


where 

I  is  solar  insolation, 

C  is  radiational  cooling, 

c"  =  exp[-1.5(|^I)  (tan2S)At 

p  SURF 

S  *  local  slope 

and  At  is  the  time  from  the  onset  of  either  upslope  or  downslope 
flow.  For  upslope  flow.  At  is  the  time  from  sunrise  and  for 
downslope  flow  it  is  the  time  from  two  hours  before  sunset.  The 
factor  291.5/Tsurf  accounts  for  the  actual  slope  temperature 
which  was  a  constant  291. 5°K  in  Ryan's  case  studies. 

The  value  of  radiational  cooling,  C^,  according  to 

Ryan  is 


C  =  .75*Pa*(l  -  %cloud  cover).  (40) 

n  u 

where  the  transmissivity ,  PQ,  was  selected  to  be  0.7.  This  is 
appropriate  to  a  moist  boundary  layer  with  some  industrial  par¬ 
ticulates  (Landsberg,  1969). 

Insolation,  I,  is 
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I  *  I  *P0^S*nC*sin5*(  1  -  %  cloud  cover)  (41) 

OB 

where  IQ  is  the  solar  radiation  at  the  top  of  the  atmosphere,  C 
is  the  angle  of  the  sun  from  the  horizontal  and,  5  is  the  angle 
between  the  surface  and  the  radiation  path  from  the  sun.  More 
extensive  discussion  about  both  C  and  5  is  given  in  Ryan,  1974. 
The  above  equations  from  Ryan  are  used  to  calculate  the  upslope 
and  downslope  wind  components. 

The  east-west  component  Us  and  north-south  com¬ 
ponent  Vs  of  the  slope  wind  are, 

U  =  -W  sina,  (42) 

s  s 

V  =  -W  cosa  (43) 

s  s 

Where  a  is  the  aspect. 

Since  the  slope  wind  was  originally  formulated  to 
model  nighttime  drainage  winds  and  seems  to  do  well  in  explaining 
drainage  wind  observations,  the  downslope  wind  confidence  factor 
is  1.00.  The  confidence  in  the  daytime  upslope  wind  can  not  be 
as  high  since  the  vertical  scale  of  the  upslope  wind  is  somewhat 
greater  than  its  nighttime  counterpart.  As  a  result,  the  con¬ 
fidence  in  the  upslope  wind  is  set  to  0.75. 

2.3.3  Valley  Wind 

The  valley  or  mountain  wind  is  caused  by  tem¬ 
perature  differences  from  above  the  valley  to  above  the  plain. 
This  wind  blows  upvalley  by  day  and  downvalley  by  night.  In  the 
surface  wind  model  the  valley  wind  is  based  on  an  extensive 
series  of  observations  conducted  by  Davidson  and  Rao  (1957). 

They  concluded  that  the  valley  wind  strength  appeared  to  be  a 
function  of  valley  width  and  depth  and  blew  from  the  valley  to 
the  plain  irrespective  of  valley  slope.  Their  conclusions  are 
further  verified  by  the  observations  of  Buettner  and  Thyer  (1965) 
for  winds  in  a  mountain  valley  near  Mt.  Rainier. 
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Ryans 's  (1974)  expression  for  the  valley  wind  at  6 
meters  elevation  is 


Wv  *  (v  +  10.00  sin  ( 15t+244 )  +  1.69  sin  ( 30t+51 . 2 ) ] ‘Pg. Y  (44) 
where 

t  is  the  time  of  day, 

V  =  .89  +  sin  (.986D  +  79.88)  is  an  amplitude  that  depends 

on  day  of  year,  D, 

Y  is  average  of  left,  Y  ,  and  right,  Y  ,  slopes  to  the 

horizon.  r 

PB  is  the  transmissivity 

is  the  magnitude  of  the  valley  wind  in  the  valley  direction, 

O  .  The  east-west  and  north-south  components  are  given  by 

U  =  -W  sin  0  (45) 

v  v 


V  =  -W  cos  0 
v  v 


(46) 


This  empirical  fit  was  designed  to  yield  an 
upvalley  wind  speed  for  presumably  optimum  conditions.  These 
occur  in  a  well  defined  valley  at  the  summer  solstice.  At  mid¬ 
day  the  optimum  wind  speed  corresponds  to  observed  winds  in  a 
Southern  California  Mountain  Valley.  Equation  (44)  gave  reason¬ 
able  values  for  the  Mt.  Rainier  valley  winds  as  well.  Because 
of  the  good  agreement  that  the  model  valley  wind  had  with 
Buettner  and  Thyer's  observations,  the  valley  wind  confidence 
factor  is  1.0. 


2.3.4  Ridge  Enhancement 

Speed  up  or  ridge  enhancement  of  flow  over  hills 
and  ridges  is  a  commonly  observed  phenomenon  that  is  especially 
crucial  to  a  realistic  surface  layer  wind  model.  Jackson  and 
Hunt  (1975)  analytically  determined  the  increase  in  wind  speed 
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c  - 


for  an  adiabatic  turbulent  boundary  layer  over  a  uniformly  rough 
low  two-dimensional  hill  (hs/ls  <  .05,  where  hs  and  ls  are 
respectively  the  characteristic  height  and  length  of  the  hill.) 
Their  solution  predicts  a  "speed-up  factor"  AI  at  the  top  of  the 
hill  of 


AE  *  2hs/ls  (47) 

Bradley  (1978)  reported  observations  of  wind  flow 
over  a  large  hill  in  Australia.  The  size  of  the  hill  lay  outside 
the  range  of  validity  for  the  Jackson  and  Hunt  theory  above, 
nevertheless,  the  observed  speed-up  agreed  well  with  that  pre¬ 
dicted  by  (47).  Bouwmeester  (1978)  further  refined  the  estimates 
of  ridge  enhancement  to  include  surface  roughness  and  upwind  and 
downwind  slopes.  Hunt  (1980)  reported  several  observations  of 
speed  up  over  ridges  which  substantiates  (47).  However  for 
steeper  hills  or  ridges,  flow  separation  limits  the  speed  up  fac¬ 
tor.  Hunt  further  reported  on  the  applicability  of  extending  the 
two-dimensional  result  to  three-dimensions. 

Observations  at  Brentwood  Knoll,  a  130  m  high 
three-dimensional  round  hill  in  southwest  England,  showed  a 
speed-up  factor  of  2.2  to  2.3  as  against  a  theoretical  factor 
using  equation  (47)  of  2.0.  Other  examples  of  three-dimensional 
observations  which  agree  with  two-dimensional  results  are  given 
in  Jackson  and  Hunt  (1975). 

In  terms  of  the  terrain  wind  parameters  (TWP's)  the 
speed  up  factor  at  the  ridge  crest  is  given  as  follows:  in 
Figure  13  the  hill  height,  hs ,  is  approximated  as  twice  the  ver¬ 
tical  distance,  DH,  between  the  zero  curvature  point  and  the 
ridge  crest.  The  hill  length,  ls,  as  defined  by  Jackson  and  Hunt 
is  estimated  as  the  horizontal  distance  between  the  zero  cur¬ 
vature  point  and  ridge  crest.  In  terms  of  the  TWP's  the  speed  up 
factor  at  the  ridge  crest  is 
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(48) 


The  calculated  speed-up  is  strictly  valid  for  the  ridgecrest. 
Observations  and  modeling  of  wind  flow  over  ridges  and  hills, 
however,  show  that  faster  flow  is  not  solely  limited  to  the  very 
top  of  a  ridge  or  hill.  It  extends  over  the  entire  top  area. 

The  exact  functional  form  of  the  speed-up  factor  is  not  presently 
known.  For  purposes  of  modeling  the  variation  of  the  speed-up 
factor  from  zero  at  some  point  on  the  lower  hill  or  ridge  to  its 

full  value  at  ridge  or  hill  top,  an  exponential  function 
e-Ds/Hr 

was  constructed.  The  Ds  represents  the  distance  from  the 
calculation  point  (x,y)  to  the  ridge  or  hill  top.  The  exponen¬ 
tial  factor  models  a  decrease  in  the  speed-up  factor  from  its 
full  value  at  ridge  top  to  i/e  at  the  distance  Ds  -  Hr.  If 
Ds  exceeds  Hr,  the  speed-up  factor  is  zero. 

At  .  e-D=/HR 

hr 

where  s  =  {(  x  -  xR) 2  +  (  y  -  yR)2>1/2 


Since  a  number  of  observations  validate  the  theoret¬ 
ical  speed  up  of  ridge  enhancement,  the  nominal  confidence  in  it 
is  set  to  1.0.  For  points  away  from  the  ridge  top,  both  the 
enhancement  and  confidence  drop  off.  It  is  felt  that  an 
appropriate  expression  of  the  confidence  is  the  exponential  fac¬ 
tor  which  spreads  the  enhancement  over  the  ridge. 

2.3.5  Flow  Separation 

Separated  flow  occurs  in  the  lee  of  sufficiently 
steep  terrain  features.  Its  occurrence  and  extent  is  especially 
sensitive  to  the  steepness  of  the  feature  and  the  atmospheric 
stability.  In  addition,  separation  is  different  for  ridges  and 
isolated  hills.  In  the  case  of  flow  over  a  ridge,  the  stream¬ 
lines  form  closed  loops  and  the  mean  streamline  from  the  separa¬ 
tion  point  reattaches  at  the  downwind  end  of  the  separation 
bubble  as  shown  in  Figure  (14).  For  flow  over  a  hill,  the 
separation  is  much  more  complicated,  some  idea  of  the  complexity 
can  be  seen  in  figure  (15)  which  shows  separated  flow  over  a  hill 
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viewed  from  above  (a)  and  £rom  the  side  (b). 

Due  to  the  differences  between  two  and  three-dimensional 
flow  separation/  the  surface  wind  model  calculates  separation  for 
ridges  and  hills  differently.  The  algorithm  that  determines  if  a 
grid  point  is  within  a  separated  zone  must  first  discriminate 
whether  a  hill  or  ridge  lies  upwind  of  the  grid  point.  This 
discrimination  is  based  on  hill  TWP's.  Next,  it  must  be  deter¬ 
mined  whether  or  not  separated  flow  is  occurring,  and  finally, 
whether  or  not  the  grid  point  lies  within  the  separation  bubble. 

The  algorithm  discriminating  a  three-dimensional 
hill  from  a  two-dimensional  ridge  compares  hill  TWP's  in  three 
directions.  For  a  ridge,  the  ratios  of  hill  height  to  hill 
length,  DH/HR,  both  45°  to  the  left  and  45®  to  the  right  of  the 
grid  point  in  the  upwind  direction  must  be  at  least  75%  of  the 
upwind  value.  Otherwise,  with  a  positive  DH/HR,  a  hill  is 
assumed . 

In  the  case  of  a  ridge,  separation  is  based  on 
whether  the  boundary  layer  is  convective  or  neutral/stable.  For  a 
neutral/stable  layer  two  criteria,  for  separation,  and  full 
separation,  are  used.  The  separation  criterions: 

§2  >  0.04  +  N_2  *  500.  (49) 

HR  B 

2  2 
where  N_  is  the  Brunt-Vaisala  frequency  squared.  N  is 
B  D 

based  on  the  stability  and  the  surface  temperature.  The  cri¬ 
terion  for  full  separation,  or  maximum  extent  of  the  separation 
region  is 

~  >  0.06  +  N  2  *  600.  (50) 

HR  b 

The  downstream  length  of  the  separated  flow  at  full  development 
is  10  hill  heights.  With  less  than  full  development  a  proportion 
of  the  10  hill  height  distance  is  used.  This  proportion  is  equal 
to  the  ratio  of  the  DH/HR  exceedance  of  the  separation  criterion 
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to  the  difference  between  the  criteria  for  full  separation  and 
separation.  The  separation  criterion  for  a  convective  boundary 
layer  is  arbitrarily  set  to  0.05.  For  a  convective  boundary 
layer,  the  separated  flow  is  arbitrarily  considered  to  extend 
downstream  for  S  hill  heights. 

For  a  three-dimensional  hill,  separation  is 
assessed  according  to  the  results  of  Bouwmeester  (1978). 

Separation  occurs  for  hills  with  lee  slopes  greater  than  0.28  in 
all  atmospheric  stabilities.  The  extent  of  separation  is  related  to 
the  hill  Froude  number. 


Fh  a  Vnb  hr  (51> 

where  is  the  general  surface  wind  magnitude  and  HR  is  the  hill 
height  TWP.  According  to  Bouwmeester  for  greater  than  10,  and 
a  neutral/stable  boundary  layer  the  separated  flow  extends  10 
hill  heights  downstream.  Otherwise,  the  separation  distance,  SD, 
is  related  to  the  Froude  number  by 

SD  «  <3  +  Fh  *  1.21)  *  F”  .  (52) 

For  a  convective  boundary  layer,  the  surface  wind  model 
arbitrarily  sets  the  separation  distance  at  five  hill  heights. 

When  a  point  is  found  within  a  separated  flow 
region  based  on  the  hill  TWP’s  upwind,  a  diagnostic  message  is 
ouput  and  the  confidence  factor  in  the  calculated  wind  reduced. 
This  reduction  is  effected  by  adding  the  general  wind  magnitude 
to  the  denominator  in  the  overall  confidence  factor  calculation 
(see  next  section).  In  all  cases,  the  relative  confidence  factor 
in  separated  flow  is  zero,  but  no  change  in  the  calculated  wind 
is  made. 

2.3.6  Resultant  Wind  and  Confidence  Factor 

For  each  grid  point  the  relevant  wind  perturbation 
for  sheltering  and  diverting,  slope  wind,  valley  wind,  ridge 
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enhancement,  and  separation  are  calculated  and  decomposed  into  U 
and  V  components.  At  this  time  the  perturbation  due  to  separa¬ 
tion  is  ignored.  The  magnitude  and  direction  of  the  resultant 
wind  after  accounting  for  the  other  four  components  is  calculated 
as 


w  = 


4 

I 

i  =  l 


z  V. 
i  =  l  1 


1/2 


and 


(53) 


tan 


-1 


z  u. 

1 


fv‘ 


(54) 


=  surface  wind  components  after  accounting 
for  sheltering  and  diverting 
»  perturbation  components  due  to  slope  wind 
»  perturbation  components  due  to  valley  wind 


where 

VV2 

U3.V3 

U4,V4  =  perturbation  components  due  to  ridge  enhancement 

A  resultant  confidence  factor  is  also  derived;  it 
includes  an  additional  component  assigned  in  case  of  flow 
separation.  The  expression  for  the  confidence  factor  is 


Confidence 


«  Z  CF.SP.  / 
i-1  1  1 


5 

Z  SP. 

i-1  1 


(55) 


where 


CP1 

s 

confidence 

in 

CF2 

at 

confidence 

in 

cf3 

- 

confidence 

in 

CF4 

at 

confidence 

in 

cf5 

m 

confidence 

in 
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4 


and 

SP^  »  Speed  of  sheltering  -  diverting  wind 
SP2  a  Speed  of  slope  wind 
SPj  *  Speed  of  valley  wind 
SP^  =  Ridge  enhancement  speed 

SP_  =  General  wind  speed,  assigned  in  case  of  separation. 

1  ’ 

I 


1 
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2.4  INTERPOLATOR 


The  resultant  surface  winds  lie  on  an  optimized  nonuniform 
grid  as  derived  from  the  terrain  processing  routines.  The  use  of 
the  surface  winds  may  be  facilitated  by  conversion  to  a  uniform 
grid  whose  spacing  and  location  need  not  correspond  with  that  of 
the  nonuniform  grid.  The  objective  analysis  of  Cressraan,  (1959) 
was  selected  to  produce  a  uniform  grid  because  of  its  versatility 
and  computational  efficiency. 

Input  to  the  interpolation  routines  includes  the  location 
of  the  nonuniform  grid  points,  the  wind  components,  and  the  con¬ 
fidence  factor  at  each  grid  point.  The  geographic  location, 
extent,  and  spacing  of  the  uniform  grid  are  also  required. 
Finally,  the  radius  of  influence,  number  of  interpolation  scans, 
percent  of  the  influence  radius  to  be  used  on  each  scan,  and  con¬ 
vergence  criterion  for  an  acceptable  analysis  are  needed. 

The  objective  analysis  requires  a  first  guess  field  which 
if  good  reduces  the  number  of  objective  analysis  scans  necessary 
for  convergence.  In  the  case  of  surface  winds,  the  first  guess 
wind  fields  for  each  wind  component  and  for  the  confidence  factor 
is  a  constant:  the  mean  of  the  nonuniforra  grid  wind  components 
and  confidence  factor  fields.  Specifically 


N 

Z  U./N 
i-1 


N 

Z  V  /N 
i-1  1 


N 

£  CF  /N  , 
i-1  1 


(56) 


(57) 


(58) 


where  N  is  the  number  of  non-uniform  grid  points. 


53 


The  first  guess  field  is  then  exposed  to  the  "observation" 
(nonuniform  grid)  points.  The  determination  of  a  uniform  grid 
value  of  surface  wind  considers  all  the  nonuniform  point  values 
within  the  radius  of  influence.  The  difference  between  a  uniform 
grid  value  and  "observation"  value  constitutes  an  interpolation 
error  C  for  each  component.  For  the  U  component,  for  example: 


uj 


Each  error  is  then  weighted  by  the  distance  between  the  uniform 
grid  point  and  "observation"  point.  The  weighting  function  for 
each  non-uniform  observation  point  is 


w. 


3 


(59) 


where  d^  is  the  distance  between  the  grid  point  and  the  obser¬ 
vation  point  and  Rf  is  the  radius  of  influence.  The  radius  of 
influence  is  generally  chosen  so  that  at  least  one  nonuniform 
point  is  within  Rf  at  any  uniform  point.  The  interpolation 
correction,  C,  applied  to  each  wind  component  (or  confidence 
factor) ,  C,  at  a  uniform  point  consists  of  the  mean  of  the  pro¬ 
ducts  of  interpolation  error  times  weighting  factor  or 


where 


C. 

3 


-W 


and  J  is  the  number  of  non-uniform  grid  points  within  the  radius 
of  influence. 
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The  correction  is  applied  to  each  component  as 

U*  »  UQ  -  C  (60) 

V'  *  V  -  C 
o 

CF'  »  CF  -  C 
o 

Applying  corrections  (60)  at  each  uniform  grid  point 
constitutes  one  interpolation  scan.  In  subsequent  scans,  the 
influence  radius  is  reduced  by  a  specified  percentage  to  enhance 
localized  features.  Between  each  scan  a  smoothing  function  is 
applied  to  each  wind  component  and  to  the  confidence  factor.  The 
smoothing  function  in  terms  of  the  U  component  is: 

U  =  1/2  U  +  1/8  Z  Ui  (61) 

where  Z  denotes  the  sum  of  the  values  at  the  four  nearest  grid 
points.  Application  of  (61)  serves  to  prevent  discontinuities 
from  arising  in  data  sparse  areas. 

Normally  a  maximum  of  five  scans  is  necessary  to  perform  an 
analysis.  After  each  scan,  the  maximum  correction  during  that 
scan  is  compared  to  a  convergence  criterion.  If  the  criterion  is 
satisfied,  the  analysis  is  terminated  and  the  results  output. 

The  steps  in  the  analysis  are  summarized  in  figure  16. 
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Now 


Local  surfaca  winds  and  confidence 
factor  (or  non-uni Cora  grid. 


First  guass  coaponants  ara  aaana  of 
input  data.  Uniform  grid  is 
darivad  by  sw  cornar  point,  tna 
aash  diaansions,  and  the  length 
and  width. 


The  subscript,  k,  refers  to  the 
uniform  grid  field. 


weighting  is  dependant  on  distance 
between  grid  points  if  witnin 
radius  of  influanca.  Influence 
radius  must  be  'tuned*  according 
to  application. 


Test  for  satisfactory  convergence. 
Convergence  occurs  when  the  aaxiaum 
correction  in  the  field  is  below 
a  tolerance  criterion.  If  not,  the 
grid  values  U^,  Vk  are  again  sub¬ 
jected  to  an  additional  correction 
scan. 


Figure  16.  Interpolation  of  Local  Surface  Winds  to  Uniform  Grid. 
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SECTION  3 
RESULTS 


This  section  presents  four  examples  of  surface  winds  calcu¬ 
lated  by  the  SWM  for  a  region  of  complex  terrain  in  Central 
Europe.  The  examples  were  chosen  to  isolate  mechanical  and  ther¬ 
mal  influences  on  the  local  wind  under  typical  weather  scenarios 
for  the  region.  In  these  examples  the  general  surface  wind  was 
specified  and  not  derived  from  the  free  atmosphere  winds. 

Examples  for  both  winter  and  summer  are  included. 

Figures  17  and  18  show  respectively  the  nonuniform  and  uni¬ 
form  output  grid  of  winds  calculated  for  a  high  value  of  the 
general  surface  wind.  For  convenience  a  template  showing  the 
terrain  in  the  example  region  is  included  with  each  set  of 
figures.  Conditions  are  for  a  winter  night  with  clear  sky  and  an 
easterly  general  wind  of  10  m/sec.  These  conditions  highlight 
ridge  enhancement  which  is  especially  evident  along  the  lower 
left  edge  of  the  figure  as  well  as  just  to  the  upper  right  of  the 
figure  center.  Sheltering  is  also  evident  in  the  lower  figure 
center.  A  careful  look  at  the  lower  right  hand  corner  will 
reveal  slight  diversion  of  flow. 

In  contrast  Figures  19  and  20  show  respectively  the  non- 
uniform  and  uniform  grid  of  winds  for  a  light  wind  condition. 

The  conditions  are  the  same  as  for  Figures  17  and  18  except  the 
general  wind  was  easterly  at  lm/sec.  With  these  conditions,  weak 
but  very  evident  slope  winds  appear.  An  especially  good  example 
is  in  the  lower  right  corner  where  slope  winds  diverge  from  a 
hill  top.  The  valley  wind  was  too  weak  in  this  case  to  make  an 
appearance.  This  is  consistent  since  a  well  defined  valley  is 
not  present  in  the  region. 

Figures  21  and  22  again  show  the  nonuniform  and  uniform 
grid  of  winds  for  a  high  wind  condition.  In  this  case  we  have  a 
summer  day  with  clear  skies  and  a  northwesterly  general  wind  at 
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Figure  17.  Surface  Wind  Field's  Uniform  Grid,  Winter 
General  Surface  Wind. 
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Figure  19.  Surface  Wind  Field,  Uniform  Grid,  Winter  Night,  Clear  Sky,  1  m/sec  Easterly  General 
Surface  Wind.  * 
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Figure  21.  Surface  Wind  Field,  Uniform  Grid,  Summar  Day,  Clear  Skies,  10  m/sec,  North  Westerly 
General  Surface  Wind. 


A  A  ■»  »  /  4  -4  *7777/7 

A  A  «*  »  ''*'*<*  A  +  *  *  ^  *  *  * 

A  a  i  *777  /  4  +  *  7  *  *  *  7 

A  A  7  7/7  h  *  //// 

A  /?  A  7/7  // 


*  A  A  *  7  7  /  7  4  A 

A  7 

<  7  /7  A 

7 

a  7  77 A 

7 

A  A  7  7  7  7  4  a 

7  7 

„A  *  7  7  7  7  a 

7  7  7 

*  777/7 

7 

7  7  /  7  / 

7  *  7  7  / 

7  7  7  7  7 

7  7*7  7 

7  7  7  7  7  7  7 


7a  777  77*7 
/  7  7  7  7 7  7 

7  77777777  7  7. 


7  7 

7  7  »  m 


7*  7777 

*  7 

* 

*  7  7  7  7  7 


*77 
i  7 
7  7 

7 

/ 


7*7 

* 


7  7 

* 

A  T  ▼  T 

» 

> 


7  7*777 


*  *  A  a  *77* 


A  7  A  A  7 

>  A  A  *  7  7*i  A 

*  r  t  ******  77777**4** 
*  A  A  7  7  7  *  a  * 

77*  *  a  7  7  a 

7  7  7 a  a  *  *  7  7  7  t  a  7  7 

/ 7  7  a  *  7  7  7  7  *  / 

7///aaaa  7777  77*777 


7 

7  7 

7  7  7  7 
7 


4  4  4  0  ~  P  4 

A  * 

A  A  i 

*  A  A  *  A  ▼ 

^  4  &  &  y  y 

A  7»y 

7  A  7  7  7  7  7 * 

7  a  7  // 


7  7/*a  * 

77  777  * 

/  /7*  AAA 

/  /?*  A  A 

7  7  7  77  7 1  a  a 

7  7  7  a 

A 

7  7  *  A  jf 

*77 7  *  *  *  * 

7  7  *  *  a 
7  7  7  7  7  /f  a  a 

7  7  *  « 

// 7  * 

7  7/7  *  A  A 

/  / /7  A 

7  7//*  « 

7  7/77*a 

77  777777 

77  7  7  7 

7  7  7  7 


7  / 

7  7 

7  7  7  7  7 

7  7  7  7 

7  7  7 

A  A  A  7 

y  t  * 


7/7  7*1  T  * 

*  '*.7*4  a« 


/  7  7  7  7 

7  7 

*  7  *  7  * 

A  7  *  * 
A  T 

A  7  V 

A  7  77  7 

*  7  7  7 7 7 


m/sec,  North  Westerly 


lOm/sec.  Figure  22  is  most  interesting.  It  very  clearly  shows 
ridge  enhancement  in  the  upper  left  corner,  just  below  middle 
along  the  left  edge,  lower  left  corner,  just  above  the  figure 
center,  and  below  and  to  the  right  of  center.  Sheltering  is  also 
clearly  evident  in  the  figure's  upper  left  corner,  lower  left 
bottom  edge,  upper  right  edge,  and  along  the  lower  right  bottom 
edge.  These  locations  compare  favorably  with  the  topographic 
features  of  the  region.  Diversion  is  best  seen  in  the  lower  and 
upper  right  corners. 

Finally,  Figures  23  and  24  are  for  light  winds  during  a 
summer  day.  The  sky  is  clear  and  the  general  wind  is  north¬ 
westerly  at  lm/sec.  The  principal  local  wind  influence  is 
upslope  winds  perhaps  best  illustrated  along  the  figures  lower 
left  edge.  Again  the  valley  wind  appears  too  weak  to  show  at 
this  scale. 
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Figure  24.  Surface  Wind  Field,  Uniform  Grid,  Summer  Day,  Clear  Skies,  1  m/sec.  North  Westerly 
General  Surface  Wind. 
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SECTION  4 
CONCLUSIONS 

The  preliminary  exercises  of  the  surface  wind  model  have 
shown  that  it  produces  generally  reasonable  and  expected  results 
with  good  computational  efficiency.  It  has  been  developed  to  a 
respectable  degree  of  credibility  as  measured  against  the 
necessary  tests  of  verification,  validation,  clarity  of  results, 
and  workability.  Some  specific  comments  with  respect  to  these 
elements  of  credibility  follow. 

4.1  VERIFICATION 

The  model  possesses  reasonable  logical  consistency.  All 
components  of  the  model  owe  their  parentage  to  parameterization 
from  tested  empirical  or  theoretical  studies.  They  have  been 
systematically  woven  together  to  yield  a  methodology  of  con¬ 
siderable  efficiency. 

4.2  VALIDATION 

The  present  model  lacks  detailed,  rigorous  testing  against 
empirical  data.  However,  the  wind  patterns,  qualitatively, 
appear  representative  of  those  reasonably  expected  over  complex 
terrain.  The  Terrain  Wind  Parameters  can  be  expected  to  be  quite 
satisfactory  based  on  tests  of  the  terrain  fitting  efficiency  of 
the  Tchebyshev  representation.  Average  fit  efficiencies  are 
typically  greater  than  90%  for  1  km  windows,  greater  than  80%  for 
2  km  windows  and  greater  than  75%  for  3  km  window  sizes. 

Standard  deviations  are  typically  a  few  percent. 

4 . 3  CLARITY 

The  output  can  be  viewed  unambiguously  and  with  a  reaso¬ 
nable  understanding  of  what  processes  are  likely  occurring.  The 
vector  plot  approach  appears  to  be  a  desirable  format  for 
presentation. 
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4.4  WORKABILITY 


The  surface  wind  model  is  cost  effective  and  can  be  used  to 
aid  in  decision  making  if  desired.  The  existence  of  an  analyti¬ 
cal  representation  of  the  terrain  surface  allows  for  efficient 
calculation  of  terrain  wind  parameters.  A  6  km  by  6  km  square 
area  using  a  1  km  window  requires  less  than  one  hour  of  CPU  time 
on  the  UDRI  DEC  VAX11/780  to  generate  TWPs  to  calculate  wind  at 
about  150  m  resolution.  Furthermore,  the  TWP  generation  is  a 
one-time-only  exercise.  Individual  wind  fields  can  then  be 
calculated  for  numerous  situations  using  less  than  five  minutes 
of  CPU  time.  This  flexibility  aids  in  its  use  as  a  decision 
making  tool. 
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